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It is notoriously difficult to apply statistical mechanics to generally covariant systems, because the 
notions of time, energy and equilibrium are seriously modified in this context. We discuss the 
conditions under which weaker versions of these notions can be defined, sufficient for statistical 
mechanics. We focus on reparametrization-invariant systems without additional gauges. The key 
is to reconstruct statistical mechanics from the ergodic theorem. We find that a suitable split of 
the system into two non interacting components is sufficient for generalizing statistical mechanics. 
While equilibrium acquires sense only when the system admits a suitable split into three weakly 
interacting components —roughly: a clock and two systems among which a generalization of energy 
is equi-partitioned. This allows the application of statistical mechanics and thermodynamics as an 
additivity condition of such generalized energy. 
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I. INTRODUCTION 

Thermodynamics is a general tool for studying macro¬ 
scopic phenomena in terms of a small number of con¬ 
cepts, such as energy, temperature, entropy, and equi¬ 
librium. Maxwell, Boltzmann, and Gibbs have related 
thermodynamics to the mechanical laws, and their statis¬ 
tical mechanics has been extended to quantum systems 
using Bose-Einstein and Fermi-Dirac distributions, and 
von Neumann entropy, which has lead to understanding 
phenomena such as black body radiation, superfluidity 
of "^He, or the stability of white dwarfs. So far, how¬ 
ever, thermodynamics and statistical mechanics have not 
found a convincing generalization in a fully general rela¬ 
tivistic context. 

The statistical mechanics of gravitationally interacting 
objects is already complicated in non relativistic approxi¬ 
mation since gravity is long-range and attractive. But the 
problem is far deeper in the full relativistic theory, when 
we take into account the fact that thermal energy can be 
ceded to the degrees of freedom of spacetime itself. The 
difficulty comes from the modifications of the notion of 
time introduced by general relativity [1]. In this context, 
no quantity concentrates all the properties that charac¬ 
terize the non relativistic notion of time on which statis¬ 
tical mechanics and thermodynamics are built. Neither 
coordinate time, nor clock times, nor proper time along 
a worldline fully behave as the time of non relativistic 
physics. As a consequence, in a general covariant con¬ 
text there is no clear and unique notion of energy, with 
all the properties of Newtonian energy, to write statisti¬ 
cal states or for thermodynamics. Another way of seeing 
the problem is to observe that if spacetime is a dynamical 
entity, a statistical state must describe also its thermal 
fluctuations — in some sense, statistical mechanics must 
take fluctuations of time into account. 

An approach to the problem based on the idea of start¬ 
ing from arbitrary statistical states and deriving notions 
of time -thermal time- and energy, from the state it¬ 
self, has been developed in references [2-6]. However, 


this line of investigation has been delayed by a concep¬ 
tual obstacle: the difficulty of characterizing the notion 
of equilibrium within this perspective. 

Here we take a different path. Standard statistical 
mechanics can be developed starting from the “ergodic 
theorem”, namely the observation that under suitable 
(“ergodic”) conditions, the time average of an observ¬ 
able is well approximated by an average over phase space. 
An equilibrium state is then defined as the phase-space 
probability measure satisfying this equality. This paper 
explores the possibility of founding covariant statistical 
mechanics on (an appropriate generalization of) this ap¬ 
proach. More precisely, we extend the concept of time- 
average to reparametrization-invariant systems by means 
of a “clock”, -a generalization of the notion of time- and 
search for conditions under which the corresponding sta¬ 
tistical state can be computed. 

We develop this idea in the context of a 
reparametrization-invariant system with a finite - 
possibly large- number of degrees of freedom and no 
other gauges beside coordinate-time reparametrization 
invariance. This is sufficient to show that a truly 
Newtonian time, with all the properties of time in the 
non relativistic theory, is not necessary to use statistical 
mechanical methods. Generalization to field theory and 
systems with additional gauges is considered elsewhere. 

We report two main results. First, when the system 
admits a suitable split into two subsystems, the split itself 
allows defining a generalized time average and a statisti¬ 
cal state. Roughly speaking, one subsystem can be taken 
as a “clock system” to describe thermal properties of the 
other. Statistical mechanics relies on this split, not on a 
peculiar choice of time variable in the clock system. 

The structure defined above, however, is not sufficient 
to make sense of equilibrium. Our second result is that 
a split of the system into three subsystems is enough for 
having a notion of equilibrium. Roughly speaking: a 
clock-system and two systems in equilibrium with respect 
to each other. If the split satisfies certain conditions, the 
full machinery of equilibrium statistical mechanics can 
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be used, and there will be an equilibrium configuration 
which maximizes entropy. The key condition is additiv¬ 
ity: the split between system and clock defines a con¬ 
served quantity which is the sum of one component per 
subsystem. 

The paper is organized as follows. Section II reviews a 
formulation of Hamiltonian mechanics sufficiently wide to 
include general covariant systems, following Dirac [7, 8 ], 
and Souriau [9]. We call this generalized Hamiltonian 
formalism “relativistic”, following [10]. Section HI intro¬ 
duces the notions of time average, statistical state, and 
ergodicity in this relativistic framework and defines the 
microcanonical ensemble when the system splits into two 
components. Section IV derives the notion of equilib¬ 
rium between two systems with respect to a third one, 
which yields to the definition of entropy, and suitable gen¬ 
eralizations of the notions of energy, temperature, and 
time. In the Appendix C we illustrate some aspects of 
the formalism by applying it to a radiation-dominated 
Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmol¬ 
ogy- 



Figure 1. Generalized hamiltonian mechanics 


Equivalently, Y = Yc\t,, where Yq is the vector field, on 
V, satisfying the (generalized) Hamilton equations 


II. RELATIVISTIC MECHANICS 

Classical mechanics describes how observables evolve 
in Newtonian time. This representation of dynamical 
phenomena is not possible for general relativistic sys¬ 
tems. In this case, mechanics describes correlations be¬ 
tween partial observables (which include time, or “clock” 
observables), i.e., it describes how partial observables 
change with respect to one another. 


i^x{Yc) = -dC. (2) 

Different constraints vanishing on the same surface define 
locally rescaled Y and therefore different parametriza- 
tions of the same orbits. Equivalently, Y is deter¬ 
mined by ( 1 ) up to a local rescaling, which affects the 
parametrization of the orbits, and has no physical mean¬ 
ing. A state defines correlations between the partial ob¬ 
servables (functions from E to K). The physics is fully 
coded by the presymplectic space (E,u; 5 ]). 


A. Dirac’s generalized Hamiltonian mechanics 

Systems like a free special-relativistic particle or 
the gravitational field in general relativity can be de¬ 
scribed by a manifestly covariant action, invariant under 
reparametrization of the evolution parameter. A Legen¬ 
dre transformation leads to the symplectic space {X, ujx), 
called the relativistic, or extended, phase space, and a list 
of vanishing constraints {Ci = 0 ; i = l..k}. The con¬ 
straints define a presymplectic surface (E,a;s = wxls): 
of codimension k. The dynamics of the system is fully en¬ 
coded in the fc-dimensional orbits in S, called states (or 
motions), which are the integral surfaces of the null di¬ 
rections of ujs- The set F of these states is called physical 
phase space and is naturally equipped with a symplectic 
structure defined by ws = TT*u}r, where E A F (see figure 
1). For a detailed review of this presymplectic formalism, 
see [ 10 , 11 ]. 

In what follows, we restrict to systems with one¬ 
dimensional orbits (fc = 1 ), i.e., no gauge except 
reparametrization invariance. The orbits are the integral 
curves 7 of a vector field Y, on E, satisfying 

( 1 ) 


B. Non relativistic examples 

Conventional Hamiltonian mechanics can be reformu¬ 
lated in this “relativistic” language. A non relativistic 
system defined by a non relativistic phase space (FqjWo), 
a non relativistic Hamiltonian Hq, and Newtonian time 
t is equivalent to a relativistic system defined by 

E = Kxro , uJs = —dHo A dt + ojQ- 

By introducing a new variable pt, as the conjugate mo¬ 
mentum of t, one can also build an extended phase space 

X = T*M X Fq , LOx = dpt A dt + ujQ, 

and, with the constraint C = pt + Hq = 0, the equations 
of motion (2) reduce to usual Hamilton equations. The 
physical phase space F is isomorphic to Fq, so 

E ~ K X F. (3) 

This is called deparametrization and allows dynamics to 
be interpreted as motion in phase space (Schrddinger’s 
picture; see figure 2). 

There are relativistic dynamical systems that can be 
“deparametrized”, namely formulated also in the con¬ 
ventional non relativistic formalism. A free particle in 


a;E(y) =0. 
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Figure 2. Structure of the presymplectic surface of any non- 
relativistic system. 


Minkowski space M is a typical example: X = T*M and 
C = + w? = 0. Deparametrization requires breaking 

Lorentz symmetry, picking one Lorentz frame, where the 
constraint becomes pt + \/p^ + rn? = 0, and a Lorentz 
time t = as a special independent evolution parame¬ 
ter, even if physically there is nothing special about it. 

But there are also dynamical systems that do not ad¬ 
mit a deparametrized version, namely that cannot be ex¬ 
pressed in the conventional Hamiltonian language. An 
example is given below. 


C. A relativistic example 

Consider the system defined by 

X = {a, Pa, x,p} , ujx = dpa Ada + dp A dx 

with a potential V (x) constraining x inside a finite vol¬ 
ume. This systems looks similar to the classical descrip¬ 
tion of a particle in a box evolving in Newtonian time. 
However, instead of the Newtonian time term pt, there 
is the pendulum Hamiltonian ^{Pa + whose conju¬ 
gate variable, the phase </> of the pendulum, is periodic: 
tant/) = a/Pa- The system describes the relative evolu¬ 
tion of a particle in a box and a pendulum, but it cannot 
be deparametrized, because there is no globally mono¬ 
tonic time variable along the orbits: the orbits are closed. 
This system admits no conventional Hamiltonian formu¬ 
lation. 

The question we address in this paper is whether sta¬ 
tistical mechanics and thermodynamics methods can be 
extended to systems like this, that do not admit a New¬ 
tonian time parameter, or, more generally, that do not 
have a preferred time-evolution parameter. The prime 
example of such systems is general relativity. 


A. Clock and time average 


In classical mechanics, we assume measurements to be 
instantaneous, to have arbitrary precision, and to not 
disturb the system. For a non relativistic system, an 
event is like a photo of all variables and their momenta 
<f{t),Pi{t), taken at a time t. The time average of a 
quantity / is defined by 

^ 9 ^ / fid"it),P^it))dt. (4) 

def T-f+oo Z1 J_j’ 


If we treat the time variable t as a partial observable 
on the same ground as the other variables, an event is la¬ 
beled by ((f,pi,t). In relativistic Hamiltonian formalism 
it is a point on the presymplectic surface S. Repeated 
measurements correspond to a list of events, all sitting 
on the same orbit 7 . 

We now want to generalize (4) to this context, i.e., we 
want to define the average of an observable / over a list 
of events. But there is no meaning in taking averages 
along an orbit without a measure on it. 

A first crucial step is the observation that a measure 
on the orbit does not require a time variable t to exist. A 
one-form generalizing dt is required, and it suffices. Let 
us therefore call “clock” a one-form 6 * on S that defines a 
volume form along each orbit. For any quantity /:£—>■ 
M and any state 7 , we can then define the 9-average of f 
along 7 by 




def 



(5) 


In the simple example in the section HC, the phase of 
the pendulum determines a clock 9 = which 

is well defined on the orbit even if 9 = d(j) makes sense 
only locally, and there is no global time parameter along 
the orbit. The operational procedure for measuring the 
average (5) is well defined. 


B. Statistical state 

Given a clock 9 and a state 7 , 0-averaging along 7 is a 
linear functional 

/ fid, a) g R- 

We call 9-statistical state the (probability) measure pg^-y 
over naturally defined by 

/ fdp0,y = fi9,-f). (6) 

Js def 


^ If one restricts to functions / : A ^ R only (where, for example, 
A C S or there is a projection E ^ A) then the statistical state 
is a measure on A. 
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So far the statistical state is nothing more than an¬ 
other name for 0-average along 7 . But in the next para¬ 
graph, we show how this statistical state can be com¬ 
puted without solving the dynamics, if the system is, in 
some appropriate sense, ergodic. 


C. Ergodicity 


Ergodicity can be understood in various manners^. A 
general framework for expressing it is the following: a 
measured dynamical system {V,Ut,vY is ergodic if for 
almost any Pq £V, 

f f{UtPo)dt= [ fdv. (7) 
T ^+00 11 J Jy 

It is common to say that statistical mechanics is valid 
when the “ergodic hypothesis” is satisfied. Concretely, 
ergodicity is notoriously hard to prove mathematically 
and notoriously very common in nature. Here we are 
not interested in proving ergodicity in a general covariant 
context, but in understanding how it can be formulated 
in such a context. 

Intuitively, (7) means that (almost) any history ex¬ 
plores V entirely and therefore the time average does 
not depend on the initial conditions Pq. Then there is 
a natural way of extending the notion of ergodicity to 
the relativistic hamiltonian context: First, we say that 
a partial observable / on E is 9-ergodic if its 0 -average 
does not depend on the state 7 of the system, i.e., 

= m Vyer. (8) 


Then, a system is called 0-ergodic if every partial observ¬ 
able is 0-ergodic. Equivalently, if the 0-statistical state 
is independent on the (pure) state 7 of the system, i.e., 


= Me V 7 e r. (9) 


In Appendix A, we prove that if the relativistic system 
(X,uJx,C) is 0-ergodic, then the statistical state pLg on 
E is given by the explicit formula 


^ e{Yc)5{C)diix 
/ e{Yc) 5{C) dpx ’ 


( 10 ) 


where Yq is defined in ( 2 ). 

An essential point here is that in spite of the appear¬ 
ance of C in this expression, this statistical state does 
not depend on the choice for the constraint C: it is en¬ 
tirely determined by the presymplectic structure. In this 
sense, it is invariant under the reparametrization gauge. 
On the other hand, it depends on the choice of the clock. 
In the next section we see how this choice can be physi¬ 
cally motivated. 


^ As in the classical works of Maxwell [12], Boltzmann [13], von 
Neumann [14, 15], Birkhoff [16] and Khinchin [17]. 

® {V, u) is a probability space and f/t : E —i E is a measure¬ 
preserving time evolution operator. 


D. Microcanonical ensemble 

Thus far, the clock has been taken to be arbitrarily 
chosen. Different clocks lead to different statistical me¬ 
chanics. We now focus on a situation when this arbi¬ 
trariness is reduced and the choice of the clock has a 
clear physical ground. 

The key requirement we demand to a realistic physical 
clock is to interact as little as possible with the system we 
measure. If we are interested in the average temperature 
on Earth, we better not have a clock that shuts itself off 
during summer. In a general covariant context, this re¬ 
quirement is captured by demanding that the global sys¬ 
tem splits into two non interacting parts—one of which 
playing the role of “clock system”—in the following sense. 

A mechanical system S, given by (EjWs), splits into 
two non-interacting subsystems 3°" and S^, if it can be 
written as the presymplectic surface in an extended phase 
space {X = A“ x X^,u!x = wx" + wxi>), defined by a 
constraint of the form 

C = C‘^ + C^ = 0, 

where C°‘ : X°‘ ^ R {a = a, b). 

Such a split determines a foliation of the presymplectic 
surface 

E= □ (11) 

where X°‘ = |Jj„ E“c and /“ is the value of (7“ on the 
leave. Each (E“o,a;x“|E“„) is the presymplectic space 
that would describe the subsystem 5“ if it was isolated, 
and we can easily check that 

dimT = dimTja -l- dimT^b -|- 2. 

The degrees of freedom of S are those of 5“, plus those 
of S^, plus one degree of freedom at the boundary of the 
two subsystems. 

Example. According to IIB, any non-relativistic sys¬ 
tem with N degrees of freedom splits into a time part 

X°‘ = {t,pt} , a;“ = dpt A dt , 

C‘^ = Pt = 

def 

and the non relativistic phase space 

= U\Pi} , w'’ = dp, A dq" , 

C\= Hoiq\p,) = I\ 

def 

where I = G K is the standard energy. 

Notice that the time part alone has zero degrees of 
freedom, and it admits a single orbit: the absolute flow 
of Newtonian time, which appears as a gauge. The other 
part has iV — 1 degrees of freedom: the dynamics at fixed 
energy, again with evolution as a gauge. By coupling the 
two systems, a novel physical degree of freedom arises, 
describing the relative evolution of the two gauge pa¬ 
rameters [18], i.e., the coupling introduces the evolution 
of the system with respect to the Newtonian time t. 
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Example. The relativistic system IIC splits into a 
pendulum component 


X°- = {a,Pa} , =dpaAda, 

C\=Upl + a^)=I‘^ 

def A 


and a particle component 


Cb 


= {a?,p} , = dp A dx, 




= ^ + V{x) = 

def 2.UI 


Once again, the pendulum part (X“,C“) alone defines 
a system with no degrees of freedom, as it has a single 
orbit. The coupled system does not describe the par¬ 
tial observables of the particle as functions of absolute 
Newtonian time, but instead their correlations with pen¬ 
dulum variables. 

Now we return to statistical mechanics. The Newto¬ 
nian time subsystem in the first example and the pen¬ 
dulum in the second are special cases of what we call a 
“clock subsystem”. In general, however, a clock subsys¬ 
tem does not need to be one-dimensional. 

Say we are interested in measuring partial observables 
of the subsystem only (i.e., f : X^ ^ M), using a 
clock in the subsystem 5“ (i.e., d is a 1-form over X“). 
Then, assuming 0-ergodicity, we get a family of statistical 
states, labeled by the generalized energy / = /*' = — J“, 
defined on Ej. In Appendix B, we prove that this state 
is given explicitly by 


^ 6iC»-I)dpx^ 

JS{Cb-I)dpx^' 


( 12 ) 


These states form the microcanonical ensemble for the 
subsystem 5^. 

Remarkably, these states depends only on the splitting 
of E. They do not depend on the specific choice of clock 
one-form on 5“, neither on the form of the constraint. 

This is our first main result: microcanonical states can 
be defined in a context much wider than conventional 
Hamiltonian mechanics. They can be defined for any rel¬ 
ativistic Hamiltonian system admitting a split into two 
non interacting components, and they depend only on 
such a split. Assuming ergodicity they allow to express 
averages along the orbits (which require solving the dy¬ 
namics) in terms of phase-space averages (which do not 
require solving the dynamics). 


IV. GENERALIZED THERMODYNAMICS 

The three laws of thermodynamics are all in question, 
in a general covariant context. Energy conservation re¬ 
quires energy to be defined, but in a generic general co¬ 
variant system there is no non trivial conserved energy: 
a covariant generator of evolution along the orbit is the 
constraint, but it vanishes identically. The second law is 


tied to a notion of equilibrium and time variable, which, 
as we see below, requires more structure to be defined. 
Finally, the “zero-th” law, namely the equality of temper¬ 
ature at equilibrium appears to be violated already for 
systems in a fixed relativistic background, as with the 
Tolman-Ehrenfest effect [19]. Under which conditions 
do we recover thermodynamics in a reparametrization- 
invariant system? 


A. First law 

We have seen that a quantity I, conserved along the 
evolution, exists if the covariant system splits into two 
non interacting subsystems. As we saw in the example 
above, I can be identified with the conventional energy of 
the non relativistic formulation, if this exists. Therefore 
/ is a generalization of the notion of energy. Given a split, 
the first law generalizes immediately to / conservation. 


B. Zero’th law 

Consider two pieces of metal at different temperatures. 
Put them into contact and isolate the full system. Ther¬ 
modynamics predicts that energy will flow between the 
two, until equilibrium is reached. Entropy measures how 
well shared is the total energy between the two pieces. 
How do we recover this prediction in the relativistic 
Hamiltonian context? 

Notice that time does not play a significant role. What 
matters is the existence of a quantity (energy) whose to¬ 
tal amount is fixed and which can be shared between sub¬ 
systems. In section HI we showed that a split of a system 
into two non interacting components was sufficient to de¬ 
fine a statistical state for one component, but a split into 
three (here the two pieces of metal and Newtonian time) 
is required in order to talk about equilibrium between 
subsystems. 

For a general covariant system, assume that S splits 
into two subsystems 5“ and Choose a clock in 5“ 
and assume also that is in turn comprised of two 
interacting subsystems and i.e., 

A = A“ X a'’ X A^ 

iox = T 

c = C‘^+ C^+ C‘^+ V^‘= = 0 , 

E = y E?„ X EjL, 

where depends on both A^ and X'^ but not on the 
clock subsystem A“, and represents an interaction terms 
that allows (a suitable generalization of) energy to flow 
between the two subsystems and 5^. 

Assuming ergodicity for the isolated system the 

microcanonical state (12) is given by 

^ S{C'^ + C- + V»--I»^)dfix^^xp 

f 5{C^ + C<^ + dfxx^xxp' 
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We can now apply the usual machinery of statistical 
mechanics to define equilibrium between the two systems 
and 5'^, i.e., if is a small perturbation (Boltz¬ 
mann’s thermalizing “grain of sand”) then 


and the microcanonical state factorizes 

cx ISiC’^ - I^)d^ix»] [<5(C^ - (/'''= - 

ibc 

and are not conserved, but they sum up to the 
total conserved I-energy quantity I = +1‘^. The 

probability distribution of having /-energy partitioned as 
(/^, /“), is given by 

n{i\p) cx n^r), 

where f2“(/“) = J 5{C°‘ — I°‘)dfj,x‘> (a = h,c) is the 
“phase space volume” of the /-energy surface. 

We can then define, for each subsystem, its entropy 


s{i) = kB\ogn{i) 

def 


and its I-temperature 


^ _ dS{T) 

Tj def dl 


(13) 


(14) 


The equilibrium value of {1^,1^) is obtained by maxi¬ 
mizing the entropy of the isolated system 




or, equivalently, by equal /-temperatures 

rpa I _ rph I 

^ I \eq — I \eq- 

The 0*^ law of thermodynamics is recovered in term of 
/-temperatures. 

A conventional hamiltonian system written in the rel¬ 
ativistic formalism (section IIB) splits into a time part 
and the system itself. For such a system, (12) reduces to 
the usual microcanonical state of non relativistic system. 
If the system splits into two weakly interacting subsys¬ 
tems, i.e., Ho Hq + Ho, then (13) and (14) are the 
usual entropy and temperature. 

But in general the /-temperature is not necessary what 
would be measured by a standard thermometer. A good 
example of this discrepancy is the Tolman-Ehrenfest ef¬ 
fect: the temperature T measured by a standard ther¬ 
mometer in an equilibrium configuration of matter on a 
fixed static gravitational field is not constant in space. 
The reason is that T is the derivative of the entropy with 
respect to the local matter energy, which is the generator 
of proper-time translations, and is not conserved, because 
of the interaction with gravity. A conserved quantity I 
is given in this case by the generator of the global time 
translations in the static gravitational field. The ratio 
between T and Tj is different in different space points, 
because the ratio between the global time and proper 
time is affected by the local gravitational field. The gen¬ 
eral formalism developed here brings back the Tolman- 
Ehrenfest’s apparent violation of the zero-th principle 
into standard (but generalized) thermodynamics. 


C. Second law 

Here we are only concerned with the problem of ex¬ 
tending the second law to the relativistic case, and not 
to address the numerous subtleties and questions raised 
by the law. 

The physics of the second law is captured by the fol¬ 
lowing fact: for most states that at time to have an en¬ 
tropy S(to) lower than the maximal entropy, dynamics 
takes the system to a state with higher entropy at time 
to + dt. (For most states, S{to — dt) is larger than S{to) as 
well.) This behavior is obtained as soon as the system is 
ergodic and the probability distributions of macroscopic 
observables are very peaked. 

Remarkably, the width of probability distributions is 
simply related to the notion of subsystem: if the system 
is macroscopic (i.e., can be split into a very large number 
of weakly interacting subsystems) then we can expect the 
second law to hold. This is a physical consequence of the 
law of large numbers and has nothing to do with the 
relativistic/non relativistic aspects of the system. 

Thus, the second law generalizes immediately to 
the relativist context: entropy increases along the 
evolution for most low-entropy configurations of a 
reparametrization-invariant system formed by a large 
number of weakly interacting components. 

D. Gauge invariance and additivity 

We have seen that the existence of a split of a 
reparametrization-invariant system into two non inter¬ 
acting subsystems implies the existence of a conserved 
quantity I. It is important to note, however, that this 
quantity is defined by the split only up to reparametriza- 
tion /—)•/ = /(/). 

Recall, indeed, that / is defined as the value of 
when the total constraint has the form = —(7“. But 
the same dynamics can be obtained from the constraint 
f{C^) = for a rather generic (bijective) function 

/. Therefore I = f{I), which is of course a conserved 
quantity as well, can be equally well taken as a general¬ 
ization of the notion of energy. 

Accordingly, we can use to define an evolution pa¬ 
rameter r along the orbits by dA/dr = {A, C^}, but this 
choice is gauge dependent. The split between clock and 
the rest is sufficient to define the microcanonical states, 
but is not sufficient to determine a gauge-invariant pre¬ 
ferred generalized time variable. In other words, the sta¬ 
tistical mechanics we defined in section III does not re¬ 
quire a preferred notion of time. 

The situation changes when we consider a split into 
three subsystems and need all systems are weakly inter¬ 
acting. The reason is that if we now replace the gen¬ 
eralized energy / = with / = /(/), its additivity, 
I = + H, is lost, unless / is affine. Therefore the I- 

energy is now defined only up to affine transformations: 
changes of units and energy zero point. 
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The crucial property for the definition of equilibrium, 
therefore, is the additivity of the conserved quantity that 
exists if the system admits a three-way split. 

If we now use the generalized energy I = to 
define an evolution parameter r along the orbits by 
dA/dr = {A,C^^}, the resulting evolution parameter r 
is no longer gauge dependent: it is the unique parameter 
(up to rescaling) whose conjugate energy is additive. (A 
moment of reflection shows that r, defined in this man¬ 
ner, is precisely the thermal time of the equilibrium state 
[ 2 - 6 ].) 

V. CONCLUSION 

We studied the conditions under which concepts and 
tools of statistical mechanics and thermodynamics ex¬ 
tend to reparametrization-invariant systems. 

To this end, we studied under which conditions suitable 
generalizations of the notions of time, energy, and equi¬ 
librium can be found, sufficient for statistical mechanics 
and thermodynamics to become meaningful. 

The key idea was to study what is needed to extend 
the notion of ergodicity (i.e., suitable independence of 
averages from initial conditions) to a reparametrization- 
invariant system. 

We introduced the following main ideas: 

Clock: Non relativistic time sets the dt measure needed 
for time averages. In a general covariant context, 
the clock describes a physical devise measuring av¬ 
erages along the evolution. 

Split into two subsystems: A split of the system into two 
non interacting parts is sufficient for determining 
a class of equivalent clocks: one-forms depending 
only on the clock subsystem define the same statis¬ 
tical state of the rest of the system. The split into 
two non interacting parts also implies the existence 
of a conserved quantity, denoted /-energy. 

Split into three weakly interacting parts: If the system 
splits into three (or more) weakly interacting parts, 
equilibrium can be defined by maximizing the en¬ 
tropy under /-energy exchange, /-energy additivity 
picks a preferred evolution parameter that can be 
identified with the thermal time of the equilibrium 
state. 

I-temperature: /-temperature (in general distinct from 
conventional temperature) is the (inverse) deriva¬ 
tive of entropy by /-energy and is constant at equi¬ 
librium. 

We showed that a split of the system into two non 
interacting parts is sufficient to define the microcanoni- 
cal ensemble and therefore to use the tools of statistical 
mechanics. 

A split of the system into three (or more) weakly in¬ 
teracting parts is needed for useful notions of equilibrium 


and temperature, and is sufficient for deriving a general¬ 
ization of all three laws of thermodynamics. 

In Appendix C we illustrate the formalism in the con¬ 
text of the reparametrization-invariant system describing 
a homogeneous FRLW cosmology filled with non homo¬ 
geneous electromagnetic radiation. We show that the 
equilibrium state is characterized by an /-temperature 
Tf constant in time. The conventional temperature T is 
related to T/ by a time-dependent scaling factor, giving 
the right time dependence of the radiation-filled FRLW 
cosmology. 

The formalism introduced in the paper tackles the con¬ 
ceptual issue of understanding thermodynamical proper¬ 
ties in the absence of Newtonian time and allows a eo- 
variant treatment of simple relativistic systems, such as 
the cosmological example in Appendix C. Is it relevant 
for gravitational physics? Investigating statistical me¬ 
chanics of more involved, realistic, gravitational systems 
necessarily require more work: one should first consider 
the present approach for systems with more than one 
constraint (more gauge symmetries) and then extend it 
to a fully relativistic field-theoretic framework. 

Nevertheless, the scheme was derived to be well suited 
for a generalization to a multisymplectic or covariant 
description of reparametrization-invariant field theories: 
within this framework, mechanics of relativistic fields, on 
a 3-|-1-dimensional manifold, can be formulated in terms 
of a finite dimensional surface E, equipped with a 5-form, 
leading to 4-dimensional orbits (see e.g. [11]). The 1- 
form used to define a clock will be replaced by a 4-form, 
providing a (gauge-invariant) measure on the orbits; the 
strategy would then consist in rethinking the notion of 
ergodicity and the derivation of statistical state in this 
context. 

If this is possible, the notions introduced here could 
perhaps be used to study the statistical mechanics and 
the thermodynamics of the gravitational field. In this 
context, the relevant split could be between regions of 
spacetime, or between spacetime and some matter, or 
else. As always in physics, approximate lack of interac¬ 
tion between clock systems and the other (partial) ob¬ 
servables should be sufficient. 


Acknowledgements 

We are grateful to Alejandro Perez and Simone 
Speziale for the interesting discussions and insights. 

Appendix A: Statistical state 

In this appendix we derive the form of the statistical 
state (10). To make the proof rigorous, some additional 
mathematical hypotheses are required: the 1-form 6 to be 
with compact support, the functions / to be continuous 
with compact supports and “S{C)dp” to be understood 
as the measure defined by the disintegration theorem. 


Let T be the parameter generated by C along the 
orbits, so that Yc = d/dr. Let be coordi¬ 

nates on X such that ojx = A dq^. Then Yq = 
~ Writing 9 = a^^dq^^ + P>^dp^,, one can 

check that, for any orbit 7 , 

/ p)r< F)r< \ 

Thus, along each orbit 9{Yc) determines the ratio be¬ 
tween the clock measure and the measure determined by 
the (gauge dependent) evolution generated by C. 

8{C)dpLx defines a measure on E, invariant under the 
(unphysical) time evolution induced by C, i.e., V(/, 

I g[P{Po,T)]6[C{Po)]dpx[Po] 

= I g[Po]S[C{Po)] dpx[Po]. 

We denote 'yp„ the (unique) orbit going through a point 
Pq G S. 0-ergodicity for a function / (8) can be written 

f f9 = m [ e , vPo. 

J-ypo d'TPo 

So, by definition, 

/ (7p„)*[/0]=/(0) / (7Po)*[^] , VPo, 

Jr. Jr 

thus VPq S E, with obvious notation, 

J f[P{Po,r)]9{Yc)\p(p„,r)dT = f{9) j 6»(yc) |p(Po,r) dr. 

By integrating both sides with the preserved measure 
6[C{Po)]dfj.x[Po], we get 

At j f9{Yc)S{C)dpx=mAT J e(Yc) 6{C)dpx, 

This gives 

- ^ J f9(Yc) 6(0) dpx 

^ Je{Yc)6iC)dpx ’ 

from which (10) follows immediately. The d{Yc) factors 
compensates the scaling of the delta function under a 
change of C, so that the state is gauge invariant. 


Appendix B: Microcanonical ensemble 

In this section we derive the microcanonical ensemble 
(12). Let’s write 

C' = C'“-fC'’ = 0,E= y 

/»+/ i >=0 


/ = = — J“ is conserved along each orbit. According 

to Birkhoff’s ergodic theorem, it is known that ergodic- 
ity can be satisfied, at most, on slice of constant I, so 
0 -ergodicity for a function f : Y R should be under¬ 

stood as 

/(0,7) = no, I) , V 7 C E, = E“, X 
Then, VPq S E/, 

J f[P{Po,r)]e{Yc)\p(p„,r) dT = f{0,I) J 6 »(yc)|p(Po.r) dr. 

By integrating both sides with the preserved measure 
S[C^{Po) - /]i5[C'(Po)]dpx[Po], we get 

f(0 n = U0{Yc)5{C^-mC)dpx 
^ Je{Yc)nC^-I)S(C)dpx ' 

The split leads to 

dpx = dpx-dpxb, 

S{C) = SiC‘^ + C^), 

Yc = Ycc+Yc.. 

As the clock 9 depends only on the partial observables of 
the subsystem 5“, diY^b) = 0. In addition, the function 
/ depends only on the partial observables of the subsys¬ 
tem so the integral 

J f9{Yc)6{C^-I)6iC)dpx 

factorizes into 

(^J 9{Yc^)S{C- + I)dpxJj (^J f6iC^-I)dpx»y 

Finally, taking the normalization factor into account, 

^ J f6{C» - I)dpx» 

’ J 6{C>> - I)dpx» ’ 

which implies immediately ( 12 ). 


Appendix C: An application: electromagnetic field 
in an expanding universe 

We illustrate some of the notions that we have intro¬ 
duced by studying a simple infinite-dimensional relativis¬ 
tic system: a FLRW universe filled with a non homoge¬ 
neous electromagnetic field (cosmic radiation) [20]. The 
dynamics is not given by the full Einstein equations since 
only the scale factor a{t) interacts with the electromag¬ 
netic field, the other degrees of freedom of the gravita¬ 
tional field are assumed to be frozen in the homogeneous 
configuration. The aim is to describe the statistics of the 
cosmic radiation. 
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The system is described by the single constraint 


2jtGpI 

3F a 


Wk 


a + 7J = 0 


and the symplectic form 


Fourier transform, the constraint can be written as 

2nGpl 




k,E{k),A{k) 


3V a a 
If we multiply it by a the constraint splits 
2ttG 


d^k = 0 . 


ay 


■Pa 


H 


k,E{k),A{k) 


d^k = 0 . 


Lu = dpa A da + 


(^dE\x) A dAi{x)') d^x, 


where 


The radiation part, in turn, splits into a sum over modes. 
Assuming a weak perturbation that allows interaction 
between the modes, the system matches the framework of 


section IV, with I^ = El 


k,E{k),A{k) = = 

—la fixed. Then, neglecting quantum aspects, one gets 
equipartition for I, i.e., Rayleigh-Jeans distribution 


ds^ = —dt^ -\- cjijdx^dx\ 

^ = J 

H = j TQo[g,E,A]^d^x = ^ j H{x)d^x, 

H{x) = g-^/^~g,,E^E^ + 

Assuming flat space (k = 0, gij = Sij) and performing 


dll, = SnkBTiu^dv. 

where Tj is a conserved generalized temperature. The 
split of the constraint generates evolution in thermal 
time. This is scaled with respect to proper time because 
proper-time evolution is generated by the conventional- 
energy H — The relation between the time- 

independent temperature Tj and the conventional tem¬ 
perature is T = which gives the time dependence of 
the temperature of the radiation dominated FLRW uni¬ 
verse. 
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